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ABSTRACT 



Context. The 2004 December 27 giant flare from SGR 1806-20 produced a radio nebula that was detectable for weeks. It was observed 
at a wide range of radio frequencies. 

Aims. To investigate the polarized signal from the radio nebula at low frequencies and to perform precise total intensity measurements. 
Methods. We made a total of 19 WSRT observations. Most of these were performed quasi simultaneously at either two or three 
frequencies, starting 2005 January 4 and ending 2005 January 29. We reobserved the field in 2005 April/May, which facilitated an 
accurate subtraction of background sources. 

Results. At 350 MHz, we find that the total intensity of the source is lower than expected from the GMRT 240 MHz and 610 MHz 
measurements and inconsistent with spectral indices published previously. Our 850 MHz flux densities, however, are consistent with 
earlier results. There is no compelling evidence for significant depolarization at any frequency. We do, however, find that polarization 
angles differ substantially from those at higher frequencies. 

Conclusions. Low frequency polarimetry and total intensity measurements provide a number of clues with regard to substructure 
in the radio nebula associated with SGR 1806-20. In general, for a more complete understanding of similar events, low frequency 
observations can provide new insights into the physics of the radio source. 

Key words. SGR 1806-20 ~ low frequency radio transients 



1. Introduction 

The 2004 December 27 flare from the Soft-y-ray Repeater SGR 
1806-20 was a major event in astronomy in a number of ways. 
■ First of all by the energy of the explosion: the brightest flash of 
radiation from beyond our solar system ever recorded. This is 
J how it caught the attention of a larger audience. Secondly, be- 
. cause the flare provided new observational data about a known 
' class of objects: magnetars, i.e., st rongly magnetized neutron 
' stars (see, e.g., iHurley et al.l l2005l) . Also, it led to specula- 
tion about a possible link with y-ray bursts (GRBs) (see, e.g., 
iTanvir et"aD [2005,) . Theorists investigated the conn ection be- 
tween the magnetic field and the explosion (see, e.g. lBlandfordI 
l2005h . Other fo cused on m odeling the fireball and th e after- 
glow (see, e.g., iNakar et al. 2005; Dai et al. 2005; Wa ng et al.1 
l2005h . Astronomers perform ed a number of follow-up observa- 
tions at various wavelengths (Rea et alj|2005t Israel et al. ||2005t 
iPalmer et al]|2005l; ISchwartz et al.ll2005l; iFender et al.ll2006l) . In 
particular, the flux from the radio nebula produced by the explo- 
sion (Gaensler et al. 2005a; Cameron et al. 2005; Taylo r et al] 
l2005h was measured very frequently in 2005 January. These ob- 
servations focused on total intensity measurements at various ra- 
dio wavelengths and on polarimetry at 8.5 GHz. Some polarime- 
try was done at lower frequencies, but without the proper correc- 
tion for the leakages (Gaensler et al. 2005b). We have performed 
accurate polarimetry at 350, 850 and 1300 MHz. Also, we were 
able to measure the Stokes I flux from the radio nebula at 350 



and 850 MHz more precisely by observing the same field again 
in 2005 April/May. In this way, we could properly subtract the 
background sources from the (u, v) data of the 2005 January ob- 
servations. We compare our measurements with those at nearby 
frequencies. 

2. Observation and data reduction 

2.1. General 

A total of 19 observations were performed in January, April and 
May of 2005. Four of these, on January 16, 20, 23 and 29 were 
alternating between 350 and 850 MHz. On January 7 and 10 
scans at 1300 MHz were also included. On January 4 we ob- 
served at 350, 650 and 1300 MHz, but the 650 MHz data was 
not used. A su mmary is show n in table [1] 

We u sed AIPS (lGreisenll2003 ) and ParselTongue jKettenis et al.1 
l2006h scripts for the reduction of all 19 datasets. The Westerbork 
Synthesis Radio Telescope (WSRT) was used for all observa- 
tions. The WSRT is a linear array with 14 equatorially mounted 
25-m dishes equipped with linear feeds. Its maximum baseline 
is 2.7 km. All datasets recorded four polarization products with 
8 IFs. 3C286 was observed before the target and 3C48 after 
RFI was excised from the spectral line data using the AIPS task 
'SPFLG'. 

Calibration was done in four steps. First we determined the vari- 
ation in system temperature as a function of time (and therefore 
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Table 1. Summary of these 19 WSRT observations 



Epoch 


Days since 


Time on 


Frequency 




burst 


source (min) 


(MHz) 


January 4 


7.6 


54 


1300 


January 4 


7.6 


91 


350 


January 5 


8.6 


462 


850 


January 7 


10.5 


107 


1300 


January 7 


10.5 


107 


350 


January 7 


10.5 


107 


850 


January 10 


13.6 


78 


1300 


January 10 


13.6 


71 


350 


January 10 


13.6 


71 


850 


January 16 


19.6 


181 


350 


January 16 


19.6 


181 


850 


January 20 


23.6 


181 


350 


January 20 


23.6 


165 


850 


January 23 


26.6 


198 


350 


January 23 


26.6 


198 


850 


January 29 


32.6 


196 


350 


January 29 


32.6 


196 


850 


April 30/May 1 


124.3 


444 


350 


May 2 


125.2 


464 


850 



also as a function of position on the sky), using the intermittent 
firing of a stable noise source. Next we performed a bandpass 
calibration using the AIPS task 'BPASS' using either 3C48 or 
3C286 or both. We applied the bandpass solution using the AIPS 
task ' SPLAT' . After that, we performed an external absolute gain 
calibration using an assumed flux for 3C48 by running the AIPS 
tasks 'SETJY' and 'CALIB'. 'SETJY' was set to use t he ab- 
solute flux density calibration determined bv lBaars et al. l (fT977h 
and the latest (epoch 1999.2) polynomial coefficients for inter- 
polating over frequency as determined at the VLA by NRAO 
staff. Finally, we self-calibrated the data for time variations in 
the relative complex gain phase and amplitude. 
Polarization calibration was performed by running the AIPS task 
'LPCAL' on 3C48 and 'CLCOR' to correct for the instrumental 
XY phase off'set. Generally, we followed the scheme for data re- 
duction of WSRT data in AIPS as outlined by Robert BraunQ, 
although we ran some AIPS tasks differently depending on fre- 
quency. Those differences mainly involved the details of polar- 
ization calibration. For instance, the leakage terms ("D terms") 
of the WSRT I Fs are channel dependent, as pointed out by 
iBrentiensI (120081 paragraph 3.2). We took account of this, by first 
averaging groups of 5 channels through the AIPS task 'SPLAT'. 
Next, we ran 'UVCOP' to make separate datasets from the av- 
eraged channels. After that, we ran 'LPCAL' and 'CLCOR' on 
each of these separately before applying the feed and XY instru- 
mental phase offset corrections by again running 'SPLAT'. 
Before imaging Stokes Q and Stokes U and before merging 
the datasets from 5 channel averaging back together through 
'DBCON', we applied a ParselTongue script for "derotation" 
to the residual data, i.e., the (u, v) data where all sources 
except the target were removed, by running the AIPS task 
'UVSUB'. The original Aips-H-i- glish script was kindly given 
to us by G. Bernardi; we modified and translated it to a 
Python/ParselTongue script on a channel by channel basis. The 
"derotation" of the visibilities is absolutely necessary, since the 
rotation measure (RM ) of SGR 1806-20 is large, 272 rad/m^ 
dGaensler et al.l |2p05a). This means that the polarized signal 
would vanish if all IFs were imaged simultaneously. For the 



350 MHz data, one really needs the derotation of the visibil- 
ities to be performed on a channel per channel basis, because 
the imaging of even one single IF would result in a severly 
corrupted measurement and underestimate of the fractional lin- 
ear polarization. The uncertainty in this RM (10 rad/m^, see 
iGaensler et alj|2005al) is too large for accurate polarization an- 
gle measurements, especially at frequencies below 1 GHz. For 
this reason we determined the RM more accurately, by fitting 
the sin 2 ■ RM ■ spectrum of either Stokes U or Stokes Q to 
its measured values at the 8 wavelengths A corresponding to the 
IFs near 350 MHz and 850 MHz. The contribution to this RM 
from the ionosphere is naturally included in this fit, at least the 
part that did not vary during the observation run. We checked 
the output of the AIPS task 'TECOR' for any significant varia- 
tions in the ionospheric Faraday rotation during every observing 
run. The ionospheric Faraday rotation computed by 'TECOR' is 
considered accurate since it does not use a model for the iono- 
sphere but actual data from the CDDIS archive. We did not apply 
the ionospheric corrections from 'TECOR' to our data because 
it implicitly assumes that one has recorded data from circular 
feeds. 

It should be clear from table[T]that the maximum observing time 
is 7.7 h due to the low declination of the source. Hence, the (u,v) 
coverage is sparse for all observations, since linear arrays like 
the WSRT ideally have 12h runs. The worst coverage was at 
three epochs when we alternated between three frequencies. 

2.2. Detailed desciption of the datasets 
2.2.1. Observations at 350 MHz 

The 350 MHz observations were performed on January 4, 7, 
10, 16, 20, 23 and 29 and April 30/May 1 of 2005. The last 
observation was made to make an accurate subtraction of back- 
ground sources possible. This mainly concerns the subtraction 
of the Lu minous Blue V a riable discussed in Supplementary 
Table 1 of IGaensler et al.l ( l2005a[) . The time resolution of all 
observations, except the first and the last was 30s. On January 4 
and April 30/May 1 the sampling times of the visibilities was 
60s. The bandwidth per IF was 10 MHz, sepai-ated 8.75 MHz 
from each other and centered on frequencies of 315.00, 323.75, 
332.50, 341.25, 350.00, 358.75, 367.50 and 376.25 MHz. The 
IFs were split into 64 channels, each 156.25 kHz wide, except 
for the April 30/May 1 observation. For that observation, the IFs 
were split into 128 channels, each 78.125 kHz wide. We used 
an automated flagger for the initial editing of our data: WSRT 
flaggefl 3C286 was included in the external gain calibration, 
along with 3C48. This was trivial, since 3C286 is unpolarized at 
this frequency. The assumed fluxes for 3C48 and 3C286 in the 
lowest frequency IF were 43.889 and 26. 106 Jy, respectively. 
The April 30/May 1 observation has the best (u,v) coverage. 
After performing 10 iterations of self calibration on this dataset 
the rms noise in the final image was 2.5mJy/beam. Its clean 
components were used to solve for the gain phases and ampli- 
tudes of the other datasets using a rather sophisticated scheme. 
First, a deconvolution of each of the 2005 January datasets 
was done in order to subtract the central region containing the 
radio nebula and the LBV, using the AIPS tasks 'IMAGR', 
'CCEDT' and 'UVSUB'. The residual data were calibrated on 
the April 30/May 1 model which had the clean components 
from the central region removed. The gain phase and amplitude 
solutions were then copied and applied to the original 2005 
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January datasets. It this way we made sure that the Stokes I flux 
from SGR 1806-20 would not be reduced by calibrating on a 
model from an observation months after the flare. As explained 
in section |2".2.3l amplitude self calibration could also reduce the 
Stokes Q flux. However, due to the large RM of the source and 
because we use 45 of the available 64 channels, the Stokes Q 
flux almost completely vanishes in a single IF at 350 MHz. Thus 
this problem does not occur, at least not before "derotation". 
PSR 1937-1-21 was observed in between SGR 1806-20 and 
3C48 for polarization calibration. This polari zation calibration 
technique is decribed in detail by iBrentiensI ^2008, paragraph 
3.2). Since the RM of this pulsar is positive. Stokes Q should 
be 90° ahead of Stokes U with increasing A^, as noted by 
iBrown & Rudnickl ( l2009l paragraph 2.3). 



Table 2. Stokes I flux measurements at 350 and 850 MHz; clean 
components from the 2005 April 30/May 1 and May 1/2 obser- 
vations were subtracted 







350 MHz 


850 MHz 


Epoch 


Days 


Stokes I 


i cr 


Stokes I 


1 cr 




since 


flux dens. 


error 


flux dens. 


error 


date) 


burst 


mJy/beam 


mJy/beam 


mJy/beam 


mJy/beam 


Jan. 4 


7.6 


186 


20 






Jan. 5 


8.6 






157 


10 


Jan. 7 


10.5 


84 


10 


97 


28 


Jan. 10 


13.6 


78 


10 


50 


22 


Jan. 16 


19.6 


16 


10 


35 


14 


Jan. 20 


23.6 


7 


10 


21 


8 


Jan. 23 


26.6 


13 


10 


17 


7 


Jan. 29 


32.6 


1 


10 


22 


9 



2.2.2. Observations at 850 MHz ("UHF high") 

We observed SGR 1806-20 on January 5, 7, 10, 16, 20, 23, 
29 and May 5 of 2005. The last observation was performed to 
make an accurate subtraction of background sources possible. 
The time resolution of all observations, except for the first and 
the last, was 30s. The sampling time of the visibilities on January 
5 and May 5 was 60s. The bandwidth of the eight IFs is 10 MHz, 
they were separated exactly 10 MHz from each other and rang- 
ing from 805 to 875 MHz. Each IF was split in 64 channels with 
a width of 156.25 kHz, except for the May 1/2 data that were 
split into 128 channels of 78. 125 kHz. The external gain calibra- 
tion was performed using an assumed flux for 3C48 of 24.240 
Jy for the lowest frequency IF. The 850 MHz were reduced in 
almost the same way as the 1300 MHz data. Only polarization 
calibration was performed slightly differently. Since the Stokes 
Q (and U) of 3C286 are not known for the "UHF high" frequen- 
cies, when the task 'CALIB' was run on this calibrator, it was 
set to solve for gain phases only and not for gain amplitudes. 



Table 3. RM measurements of SGR 1806-20. 



Epoch 


Days 


Frequency 


Measured 


1 o- 


(2005 


since 


(MHz) 


RM 


error 


date) 


burst 




(rad/m^) 


(rad/m-) 


Jan. 4 


7.6 


350 


253.24 


1.82 


Jan. 5 


8.6 


850 


253.14 


12.43 


Jan. 7 


10.5 


350 


253.65 


1.05 


Jan. 10 


13.6 


350 


261.74 


2.04 



is run on a Stokes I model, it implicitly assumes that Q=0. 
Consequently, the same model is used to derive the X gains 
from the XX visibilities as the Y gains from the YY visibilities, 
while in fact XX=I-Q and YY=I+Q, so different models should 
be used. When solving for gain phases only, the error made is 
generally considered acceptable. 



2.2.3. Observations at 1300 MHz 

We observed SGR 1806-20 at 1300 MHz on January 4,7 and 
10 of 2005. T he total intensity measu rements have already been 
published (see lGaensler et al1l2005al) . so we focused on the po- 
larized signal. However, we did check that our Stokes I fluxes 
agreed with those previously published. 

On 2005 January 4 visibilities were recorded every 60s, on 
January 7 and 10 every 30s. The eight 20-MHz IFs were cen- 
tered on frequencies of 1255, 1272, 1289, 1306, 1323, 1340 and 
1357 MHz. Each IF was split in 64 channels with a width of 
312.5 kHz. The external gain calibration was performed using 
an assumed flux for 3C48 of 17.388 Jy for the lowest frequency 
IF. 3C286 was also included in the external gain amplitude and 
phase calibration using an assumed flux of 15.550 at 1255 MHz. 
3C286 is linearly polarized. We took account of this and of the 
usual "AIPS for linear feeds" projection (R— >X,L— >Y) by plac- 
ing the assumed Stokes Q flux of 3C286 (0.594 Jy at the lowest 
frequency IF) with a minus sign at the position of Stokes V in the 
AIPS SU table. For the other IFs we kept the same ratio between 
Stokes I and Stokes Q. In this way we could use 3C286 not only 
for fixing the instrumental XY phase offset, but also for external 
gain calibration. 

Self-calibration was run to solve for the gain phases only, since 
solving for the amplitudes could reduce the Stokes Q flux. 
The AIPS task 'CALIB' cannot be set to run simultaneously 
on a Stokes I and Stokes Q model. Obviously, when 'CALIB' 



3. Results 

3. 1 . Total intensity measurements 

The total intensity flux measurements at 350 MHz were 
done by fitting a Gaussian of the same shape and size as 
the restoring beam to the (fixed) location of SGR 1806-20 
in the Stokes I images. This was d one by the AIPS task 
'IMFIT'. We used the position from iGaensler et all ( l2005al) 
(»= 18h08m39.343s,(5= -20°24'39.8") for flie fits. The re- 
sults are summarized in table |2] The error bars are conservative 
estimates from measurements of the residuals of bright sources 
in the field. The actual rms noise in these images is much lower, 
around 3mJy/beam, which is about the same as the error from 
'IMFIT'. 



3.2. Poiarimetry 
3.2.1. General 

Poiarimetry was performed on 2005 January 4, 5, 7 and 10. 
Although all of our observations recorded full Stokes, we antic- 
ipated that it would not be possible to detect the polarized signal 
from SGR 1806-20 on later dates, since the total intensity drops 
rapidly. Also, we d id not expect polari zation fractions to exceed 
the values given bv lTavloret all (120051 table 2). 
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Fig. 1. Comparison between the 240, 3 50, 610 and 1300 MH z fluxes of the radio nebula associated with SGR 1806-20. The 1300 
MHz fluxes were published previously jGaensler et al.ll2005ah . 



Table 4. Polarimetric measurements of SGR 1806-20 



Epoch 


Days since 


Frequency 


Ve^ + 


1 cr 


Polarization 


1 0- 


Polarization 


1 cr 


Reduced 


(2005) 


burst 


(MHz) 


(mjy/beam) 


error (mJy/beam) 


fraction (%) 


error (%) 


angle (°) 


error (°) 


X 


Jan. 4 


7.6 


350 


2.68 


0.81 


1.44 


0.46 


103 


39 


0.70 


Jan. 4 


7.6 


1300 


0.71 


0.66 


0.47 


0.44 


31 


26 


1.57 ' 


Jan. 5 


8.6 


850 


2.22 


0.28 


1.41 


0.20 


96 


7 


0.85 


Jan. 7 


10.5 


350 


2.30 


0.64 


2.73 


0.82 


69 


38 


0.33 


Jan. 7 


10.5 


850 


1.71 


0.52 


1.76 


0.74 


44 


9 


0.65 


Jan. 7 


10.5 


1300 


1.90 


0.37 


2.29 


0.48 


50 


7 


1.06 


Jan. 10 


13.6 


350 


1.14 


0.87 


1.46 


1.13 


36 


47 


0.79 


Jan. 10 


13.6 


1300 


1.31 


0.40 


3.65 


1.38 


69 


12 


1.28 



' Poor fit. 



3.2.2. Determining tine RIVl of SGR 1806-20 

As noted before, the rotation measure (RM) as measured by 
dGaensler et al.ll2005al 272 + 10 rad/m^) has a rather large error 
bar which translates into a polarization angle uncertainty at 1300 
MHz of 30.5°. At 850 MHz this is even 71.4°. NaturaUy, the 
RM should be determined more accurately before polarization 



angles are to be measured. 

This can be done by plotting Stokes U or Q fluxes of SGR 
1806-20 as a function of frequency and fit for the RM. We are 
in the advantageous position that these WSRT observations 
were performed with eight IFs. Over a wide span of frequencies 
there are many turns of Stokes U (or Q) since its spectrum 
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Fig. 2. Determining the rotation measure of SGR 1806-20 by fit- 
ting the sinusoidal Stokes U spectrum. Here, the fit was made 
to the values of Stokes U on 2005 January 04 at the wave- 
lengths coiTesponding to the 8 IFs near 350 MHz after the 
visibilities were "derotated" by an angle corresponding to an 
RM of -272 + 10 rad/m^. We used gnuplot to fit the function 
A ■ sin(2 ■ RM -A^ + O) for three free parameters A, RM and 0. The 
correction to the RM from this fit is -18.76 + 1.82 rad/m^. The 
reported reduced is 0.69. 
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Fig. 3. Comparison between linear polarization fractions at 350, 
850, 1300 and 8400 MHz 



o 
a 

fi 
o 
+j 

03 
N 

'P 

PL, 



160 
140 
120 
100 
80 
60 
40 
20 




OO 0.35 GHz 
O O 0.85 GHz 
□ □ 1.3 GHz 
W W 8.4 GHz 



i 



' t ♦ 



7 8 9 10 20 

Time since burst (days) 

Fig. 4. Comparison between the polarization angles at 350, 850, 
1300 and 8400 MHz 



but the latter observation yielded a very poor constraint on 
the RM, so we left it out. The 1300 MHz data also gave very 
poor constraints on the RM, thus in determining the weighted 
mean RM we ignored those, too. For the other observations, 
we plotted Stokes U per IF and solved for the RM (850 MHz) 
or the correction to the RM (350 MHz), as illustrated in figure 
|2] The results are shown in table [3] It turned out that the noise 
levels in all of the Stokes Q maps were much higher than in 
the Stokes U maps, so we did not use them. In determining the 
wei ghted mean RM we als o took into account the measurement 
by Gaensler et akl (l2005al 272 + 10 rad/m^) From this set 
of five measurements we derived an RM of 255.01 + 0.83 
rad/m^. It should be clear that, with regard to the 350 MHz RM 
measurements, the fits give the same reduced for both the 
positive and the negative correction to the initial "derotation". 
We removed those ambiguities by considering the Stokes U 
measurements near 850 MHz data on 2005 January 5. The fit to 
that data gave an RM of 253.14 ± 12.43 rad/m^ which made all 
of the positive RM solutions to the 350 MHz data very unlikely 
(^ 3.0cr level for January 4 and 7). 

It is evident that the contribution of the ionosphere to the RM, 
RMion, is included in all fits. For the 2005 January 4, 5, 7 and 10 
observations, RMio,, as reported by the AIPS task 'TECOR', is 
the range 2.1 + 0.4 rad/m^. Consequently, the interstellar RM is 
given by RMim = 255.01 - 2.1 = 252.91 ± 0.92 rad/m^. 



is sinusoidal as a function of A^. This effect is largest at low 
frequencies: at 1300 MHz, there is less than one cycle of 
A ■ sin(2 ■ RM ■ + (p), at 850 MHz there are almost two cycles 
and at 350 MHz there are 23 cycles. It is evident that the most 
accurate measurement can be made at the lowest frequency, if 
there is sufficient signal to noise. Fortunately, we could detect 
polarized signal at 350 MHz from all three observations on 2005 
January 4,7 and 10 a fter an initial "derotatio n" of our visibilities 
using the RM from iGaensler etal] (l2005al 272 rad/m^). This 
initial derotation prevents diminution of the polarized signal in a 
single IF. At 850 MHz this initial derotation was not necessary. 
The noise levels at that frequency were such that detecting a 
polarized signal was only possible on 2005 January 5 and 7, 



3.2.3. Polarization fractions and position angles 

We were able to measure the fractional Unear polarization on all 
of the four epochs mentioned in paragraph |3.2.1| At 850 MHz, 
we were not able to measure polarization on 2005 January 10. 
For the other occasions, the measured polarized fluxes, P - 
^[^^^rU^, fractions and their error bars are listed in table|4l The 
latter two quantities are depicted in figure [3] The overall conclu- 
sion is that there is no compelling evidence for any significant 
depolarization at any frequency. Only the polarization fraction 
at 1300 MHz on January 4 is low compared to the 8.4 GHz mea- 
surements, but this fraction was determined from our worst fit, 
i.e., the fit with the highest reduced 
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The polarization angles and their uncertainties are also Usted in 
table |4] The observations at 850 and 1300 MHz gave the most 
accurate position angles, with typical uncertainties of order 10°. 
They are depicted in figure|4] Here, we see compelling evidence 
for significantly different p olarization an gles with respect to the 
8.4 GHz observations from iTavlor et al.l (j2005), particularly on 
January 5 and 850 MHz and on January 10 at both 850 and 1300 
MHz. 



4. Discussion 

4. 1. Total intensity measurements 

It is clear from figure [T] that SGRl 806-20 is much dimmer at 
350 MHz than what would be expected from the GMRT obser- 
vations at 240 and 610 MHz (Cameron et al. 2005). In princi- 
ple the Luminous Blue Variable, 14" to the east of SGR 1806- 
20 (see the Supplementary Information tol Gaensler et al.ll2005al) 
should be easily distinguishable from the Soft Gamma Repeater 
in the GMRT images, even at 240 MHz. The F WHM beamsize 
reported at that frequency is 12" x 18" dChandr a 2005b). This 
makes it hard to understand the discrepancy. 
In principle the discrepancy cannot originate from the inclusion 
or exclusion of extended emiss ion. The GMRT data were cor- 
rected for this ("Chandra 2005a b). We excluded short spacings 
(< Ik/l) from our 350 MHz WSRT observations. This was actu- 
ally a necessity since these were daytime observations and solar 
interference would otherwise compromise our calibration (see 
also Brentjens 2008, end of paragraph 3.2). 
Also, it is possible that the LBV radio nebula is variable and that 
it was much brighter on 2005 April 30/May 1 than on some oc- 
casions in 2005 January. We ran the AIPS task TMFIT' on the 
map from our 2005 April 30/May 1 observation and we found a 
peak flux density of 138 + ImJy/beam and an integrated flux 
of 189 ± 2mJy/be am at the location of the LBV. The NVSS 
dCondon et al.lll99&) image of this field shows this source at the 
15 mJy level. This would indicate that the LBV has a spectral 
index of about -1.8, which is almost the index for thermal radio 
radiation. It should be noted that, at the times of the latest ob- 
servations in January 2005, when the radio nebula was relatively 
dim, there is no evidence for negative residuals in our maps that 
could be caused by the subtraction of the LBV. This indicates 
that, most likely, the LBV had the same brightness at the times 
of at least some of the 2005 January measurements as on 2005 
April 30/May 1. 

Variability at radio wavelengths of the radio nebula s from LBVs 
has been known for quite some time (see, e.g., lAbbott et alj 
11981 ). For the P Cygni nebula variability at times cales of days 
was established at cm wavelengths ( Skinner et alJ[T996) . These 
authors report a 50% increase in flux in less than two days on 
one occasion during three months of observations on every other 
day. It is unknown how these variations translate to lower fre- 
quencies. We therefore cannot completely exclude that the LBV 
was brighter at the time of the 2005 April 30/May 1 observa- 
tion than on some occasions in January 2005. Also, the spectral 
index derived above does not agree with any of the spectral in- 
dices of the four LBVs observed by Duncan & White (2002) at 
3 and 6 cm. Two of those spectral indices are close to that of 
a s pherically symmetric radially expanding stella r wind (-1-0.6, 
see lPanagia & Fellilll975l:IWright & Barlowlll975h . However, at 
these wavelengths, those systems may well be described as op- 
tically thin, which may not be the case at the frequencies we are 
considering. 

The WSRT 850 MHz Stokes I measurements are not inconsistent 



with th e 840 MHz MOST data published earlier dGaensler et al.l 
'2005a'), given the rather large noise levels in the data from both 
telescopes. The last MOST observation was taken 15 days after 
the Giant Flare (GF). Consequendy, the 850 MHz WSRT obser- 
vations after 2005 January 10 cannot be compared with other ob- 
servations in this band. The last three of the January 2005 obser- 
vations at 850 MHz were less contaminated by RFI than the first 
four, which resulted in smaller error bars on the fluxes. There is 
evidence (> 2cr level) for a deviation from a power-law decay 
from about 1 5 days after th e GF, a nalogous to the 4.8 GHz ob- 
servations bv iGelfand et al.l (I2005L paragraph 2). These authors 
also mention a gradual rebrightening from about 25 days after 
the GF, as a result of swept up ambient material. We can also see 
that in the WSRT 850 MHz data, but the evidence for this is less 
compelling, since the sampling of these observations is sparse in 
time. Consequently, it is shown only in one of our observations, 
on 2005 January 29, 32.6 days after the GF. 

4.2. Polarimetric measurements 

In figure |3] we compare the polarization frac tions as listed in ta- 
blelJwith the measurements at 8.4 GHz bv lTavloret all ( I2005L 
Table 2). In figure |4]we have done the same for the polarization 
angles. It is clear that the observations at 8.4 GHz are much more 
accurate. Still, we do not see any significant discrepancies in the 
polarization fractions. 

Our observations reveal larger polarization position angles than 
the 8.4 GHz observations. Most compelling are the observations 
on 2005 January 5 at 850 MHz and on January 7 at both 850 
and 1300 MHz. The error bar on the polarization angle at 350 
MHz on January 4 is rather large, but this measurement and the 
850 MHz measurement on January 5 show the largest differences 
with the 8.4 GHz observation, about 85°. At these times, the po- 
larization angles from the 8.4 GHz observations suggest that the 
magnetic field in the emitting plasma is aligne d preferentially 
along the axis of the radio source, on average dGaensler et al.l 
2005a). Thus, the January 4 and 5 polarization angles at 350 and 
850 MHz indicate that the magnetic field in the emitting plasma 
that causes linearly polarized radiation at these low frequencies 
is close to perpendicular to the axis of the radio source, within 
^ 20°. Possibly a different substructure in the radio nebula is be- 
ing probed. It seems hard to explain this feature without a com- 
plex model of the radio source. 



5. Conclusions 

It is striking that depolarization at low frequencies is absent. 
Also, we have shown that low frequency polarimetry of SGR 
1806-20 provides hints with respect to the detailed substructure 
of the radio nebula which cannot be derived from the extrapola- 
tion of high frequency measurements. Models for the radio neb- 
ula need to take into account a distinct source of linearly polar- 
ized low frequency radiation with magnetic fields in the emitting 
plasmas aligned quite differently from the fields that are associ- 
ated with radiation at high frequencies. 
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